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ABSTRACT 

We present the results of a deep (15 < r < 23), 20 night survey for tran- 
siting planets in the intermediate age open cluster M37 (NGC 2099) using the 
Megacam wide-field mosaic CCD camera on the 6.5m MMT. We do not detect 
any transiting planets among the ~ 1450 observed cluster members. We do, 
however, identify a ~ 1 Rj candidate planet transiting a ~ 0.8 Mq Galactic field 
star with a period of 0.77 days. The source is faint {V = 19.85 mag) and has an 
expected velocity semi-amplitude of if ~ 220 m/s (M/Mj). We conduct Monte 
Carlo transit injection and recovery simulations to calculate the 95% confidence 
upper limit on the fraction of cluster members and field stars with planets as a 
function of planetary radius and orbital period. Assuming a uniform logarith- 
mic distribution in orbital period, we find that < 1.1%, < 2.7% and < 8.3% 
of cluster members have 1.0 Rj planets within Extremely Hot Jupiter (EHJ, 
0.4 < P < 1.0 day). Very Hot Jupiter (VHJ, 1.0 < P < 3.0 day) and Hot Jupiter 
(HJ, 3.0 < P < 5.0 day) period ranges respectively. For 0.5 Rj planets the 
limits are < 3.2%, and < 21% for EHJ and VHJ period ranges, while for 0.35 Rj 
planets we can only place an upper limit of < 25% on the EHJ period range. For 
a sample of 7814 Galactic field stars, consisting primarily of FGKM dwarfs, we 
place 95% upper limits of < 0.3%, < 0.8% and < 2.7% on the fraction of stars 
with 1.0 Rj EHJ, VHJ and HJ assuming the candidate planet is not genuine. If 
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the candidate is genuine, the frequency of ~ 1.0i?j planets in the EHJ period 
range is 0.002% < Jehj < 0.5% with 95% confidence. We place limits of < 1.4%, 
< 8.8% and < 47% for 0.5 Rj planets, and a limit of < 16% on 0.3 Rj planets 
in the EHJ period range. This is the first transit survey to place limits on the 
fraction of stars with planets as small as Neptune. 

Subject headings: open clusters and associations:individual (M37) — planetary 
systems — surveys 



Introduction 



The discovery by iMayor fc Queloa (119951 ) of a planet with half the mass of Jupiter 
orbiting the solar- like star 51 Pegasi with a period of only 4.23 days shocked the astronomical 
community. The existence of such a "hot Jupiter" (H J) defied the prevailing theories of planet 
formation which had been tailored to explain the architecture of the Solar System. Since 
then, radial velocity (RV) surveys for planets orbiting nearby F, G and K main-se quence 
stars have determined that 1.2% ± 0.2% of these stars host a HJ (IMarcy et al.ll2005l . where 
a HJ is defined as a planet roughly the size of Jupiter that orbits within 0.1 AU of its star), 
with indications that this frequency depends on t he metallicity of the ho st star (such that 
the frequency is roughly proportional to lO^^^'^/'^l. Fischer fc Valentill2005l ). 



Over the last decade th ere have been nume rous surveys for extra-solar planets following 
a variety of techniques (e.g. iButler et al.ll2006l ) with the goal of determining the planet oc- 



currence rate in new regions of parameter space. A successful technique has been to conduct 
photometric searches for planets that transit their host stars. This technique is particularly 
sensitive to planets on close-in orbits. To date more than 50 planets have been discov- 
ered by this techniquj^, including numerous very hot Jupi t ers (V HJ) with orbital periods 
between 1 and 3 days. ICaudi. Seager. and Mallen-Ornelas ( 2005) used the four transiting 



plane t s discovered at the t i me by the OGLE collaboration ( Udalski et al.ll2002l : iKonacki et al. 



2003 
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Bouchv et al. 


2004 





2004 



Pont et al. 



Gould et al. 



20041 ) to determine that only 



J2006h conducted a thorough analysis 



0.1 - 0.2% of FGK stars host a VHJ. ^ , 

of the OGLE survey to determine that the frequency of VHJs is fvHj = (1/710) (1^q;54) and 



^Observations reported here were obtained at the MMT Observatory, a joint facihty of the Smithsonian 
Institution and the University of Arizona. 

^http:/ /exoplanct.eu 
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Jhj = (1/320) (lljj^), while iFressin et al.l (120071 ) found fvHj = (1/560) and fnj = (1/320). 



The SWEEPS su rvey for transiting planets in the Galactic bulge conducted with the Hubble 



Space Telescope (ISahu et al.ll2006l ) identified a putative class of ultra-short-period planets, 
or Extremely Hot Jupiters (EH J) with periods less than 1.0 day orbiting stars lighter than 
O.88M0. They conclude that ~ 0.4% of bulge stars more massive than ~ 0.44Mq are orbited 
by a Jupiter-sized planet with a period less than 4.2 days, though they estimate that this 
fraction is uncertain by a factor of 2. Note that due to their faintness the majority of the 



the T rES (e.g 



Alonso et al. 


2004 


), HAT (e.g. 


Bakos et al. 


2007 


), XO (e.g. 



20061 ) . and WASP (e.g. ICoUier-Cameron et al.l 120071 ) surveys have all discovered planets or- 



biting relatively bright stars in the Galactic field, though to date these surveys have not been 
used to calculate the planet occurrence frequency. 

While photometric surveys of Galactic field stars have been quite successful at find- 
ing transiting planets over the last few years, it is generally difficult to measure the planet 
occurrence freque ncy with these surveys (for discussions of how this can be done despit e 
the difficulties see IGouM et al.lbood : IPressin et al1l2007l : lGaudill2007l : Iseattv fc G^i^boosl ). 
The difficulty arises from the uncertainty in the parameters (mass, radius, metallicity) of 
the surveyed stars. Moreover, typical field surveys yield numerous false positives that are 
often culled in part by eye, these culling procedures are generally difficult to model in de- 
termining the detection efficiency of the survey. In contrast to field surveys, surveys of 
globular and open star clusters observe a population of stars with parameters that are rel- 
atively easy to determine en masse, moreover many of the false positive scenarios are less 
common for this type of survey. There has been significan t work inve s ted in developing 



optimum strategies to search for planets in stellar clusters ( jJaned Il996l : Ivon Braun et al. 



20051 : iPepper and Gaudill2005l ). A number of groups have completed tra r isit surveys of ope n 



clusters, including the UStAPS (Street et al. 112003 



vey by 



EXPLORE-OC (Ivon Braun et al. 1120051 ). PISCES ( Mocheiska et al. 



■ Bramich et al. 



jBurke et allbood. here after BQ6) and MONITOR jAigrain et al. Il2007h projects and a sur- 



Montalto et a. 



(iGilliland et al.l 



200G 



2005 



Hood et al. 1120051) . 
2005l . l2006h . STEPSS 



(l2007f). There hav e also been several surveys of glob ular clusters 



Weldrake et al.ll2005l : IWeldrake. Sackett fc Bridged l2008l ). 



While to date no confirmed transiting planet has been found in a stellar cluster, many 
of these surveys have placed limits on the frequency of hot transiting planets, typically as 
functions of planetary radius as well as period. The globular cluster s urveys have placed the 
most stringent constraints; the null result for the core of 47 Tucanae by IGilliland et al.l (120001 ) 
suggests that the frequency of HJ in this environment is at least an order of magnitude less 
than for the sol a r neig hborhood, while the null result for the outskirts of the same cluster by 
Weldrake et al.l (120051 ) is inconsistent with the planet frequency in the solar neighborhood 
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at the 3.3(7 level and suggests that the dearth of planets in this globular cluster is due to low 
metallicity rather than crowding effects. The open cluster surveys, on the other hand, have 
typically placed limits on the occurrence frequency that are well above the 1.2% measured by 
the RV surveys. Notably B06 conducted a thorough Monte Carlo simulation of their transit 
survey of the open cluster NGC 1245 to limit the frequency of EHJ, VHJ and HJ with radii 
of 1.5 Rj to < 1.5%, < 6.4% and < 52% respectively. The fundamental limit on the ability 
of open cluster surveys to place meaningful limits on the occurrence frequency of Jupiter- 
sized planets appears to be due to the relatively small number of stars in an open cluster. 
BOG find that for their survey strategy, ~ 7400 dwarf stars would have to be observed for at 
least a month to put a limit of less than 2% on the planet frequency, which is significantly 
larger than the typical size of an open cluster. O ne open cluster that has been a popula r 
target is NGC 6791. This cl uster is old (t 
metal rich {[M/H] hO.4, ICratton et al 



» Uy ] 
20061: 



((m — M)o ~ 12.8. IStetson. Bruntt fc Grundahll |2003| ) so that lower mai n sequence stars in 



8 Gv r.lCarraro et al.ll2006l : iKahrai et all 120071 ) 



Origlia et al. 20061 ) and contains a large 



number of stars (M > 4000M^, iKaluznv fc Udalskil ll992hT though it is also very distant 



the cluster are quite faint. It has been the tar get of three transit searches (IBruntt et al.ll2003 



Mochejska et al. II2005I : iMontalto et al.ll2007l ). the most recent of which found that their null 
result is inconsistent at the ~ 95% level with the RV HJ frequency at high metallicity. 

The paucity of stars in open clusters appe ars to limit their usefu lness as probes of the HJ 



frequency (excluding, perhaps, the result f romlMontalto et al. 



useful for probing smaller planet radii (see iPepper and Gaudi 



20071). They may, however, be 



RV sur veys have discovered a number of Neptune and super-Earth-mass p 



20061 ) . In the last several years 
anets (HN, M < 



O.lMj; 


Butler et al 


2004 


Endl et al. 


2008|: 


Fischer et a 


2007: 


Lovis et al. 


2006 


Melo et al. 


2007: Rivera et al. 


2005 




Santos et al. 


2004: 


Udrv et al. 


20061. 


2007: 


Voet et al. 


2005 


). One 


of these planets, GJ 436 b, has been discovered to transit its host star ( 


Gillon et al. 


2007). 



Little, however, is known about the frequency of these planets. Determining, or placing 
meaningful limits on this frequency would provide a powerful test of planet formation models. 



The theoretical predictions of the frequency of these o 



decline in the fre quency of HN relative to H J (llda fc Lin 



j.iects run the gamut from a steep 



Jlda fc Linll2005[ ). to HN being ubiquitous jBrunini fc Cioncoll2005l ) 



20041) , exce pt perhaps for M-dwarfs 



In this paper, the fourth and final in a series, we present the results of a survey for tran- 
siting hot planets as small as Neptune in the intermediate age open cluste r M37 (NGC 2099) 
using the M MT. We were motivated to conduct this transit survey by iPepper and Gaudi 
(I2OO5I . I2OO6I ) who suggested that it may be possible to find Neptune-sized planets transiting 
solar-like stars by sur veying an open clust er with a large telescope. The Megacam mosaic 
imager on the MMT (IMcLeod et al.ll2000l ) is an ideal facility for conducting such a survey 
due to its wide field of view and deep pixel wells that oversample the stellar point spread 
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function (PSF). Preliminary observations of NGC 6791 sug gested that finding N eptune-sized 
planets was indeed tech nically feasible using this facility (iHartman et al.ll2005l ). Using the 
formalism developed by lPepper and Gaudil (120051 ) we found that M37 is the optimum target 
for MMT/Megacam to maximize the number of stars to which we would be sensitive to 
Neptune-sized planets. We note that a drawback of this type of survey is that any identi- 
fied candidates will be quite faint making follow-up RV confirmation difficult. For planets 
significantly smaller than 1.0i?j, false positives where the transiting object is a small star 
or brown dwarf are no longer applicable. Given the depth of the survey very few giants will 
be included in the sample, and those that are, can easily be rejected based on their colors. 
Low-precision spectroscopy may be sufficient to rule out various blended binary scenarios. 
Therefore, it is reasonable to suppose that for a given small-radius candidate one could make 
a strong argument that the object is a real planet without obtaining a RV determination 
of its mass. Also note that similar difficulties will be faced by the Corot and Kepler space 
missions (albeit for even smaller planets), so our experiment may provide a useful analogy 
to these missions. 

We conducted the survey over twenty nights between December, 2005 and January, 2006, 
accumulating more than 4000 quality images of the cluster. Thi s is easily the largest te lescope 
ever utilized for such a survey. In the first paper in the series (IHartman et al.ll2008al . Paper 
I) we describe the observations and data reduction, combine photometric and spectroscopic 
data to refine estimates for the cluster fundamental parameters (t = 550 ± 30 Myr with 
overshooting, [M/H] = +0.045 ± 0.044, (m - M)y = 11.57 ± 0.13 mag and E{B - V) = 
0.227 ± 0.038 mag), and dete rmine the cluster mass function down to 0.3Mq. In the second 
paper (IHartman et al.ll2008bl . Paper II) we analyze the light curves of ~ 23000 stars observed 



by this survey to discover 1430 variable stars. In the third paper (IHartman et al.ll2008c . 
Paper III) we use the light curves to measure the rotation periods of 575 probable cluster 
members. This is the largest sample of rotation periods for a cluster older than a few hundred 
Myr, and thus provides a unique window on the late time rotation evolution of low-mass 
main sequence stars. 

In the following section we will summarize our observations and data reduction. In § [3] 
we discuss the pipeline used to remove systematic variations from light curves and identify 
candidate transiting planets. In § H] we describe the candidate transiting planets identified 
by this survey, finding no candidates that are probable cluster members. In § |5] we conduct 
Monte Carlo simulations to determine the transit detection efficiency of our survey. In §|6]we 
present our results on the limit of stars with planets for various planetary radii and orbital 
periods. Finally, we conclude in § [71 
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2. Summary of Observations and Data Reduction 



The observations and data reduction procedure were described in detail in Papers I and 
II, we provide a brief overview here. The photometric observations consist of gri photometry 
for ~ 16, 000 stars, gri photometry for stars in a field located two degrees from the primary 
M37 field and at the same Galactic latitude, an d r time-series phot ometry for ~ 23, 000 
stars, all obtained with the Megacam instrument (IMcLeod et al.ll2000l ) on the 6.5 m MMT. 
Megacam is a 24' x 24' mosaic imager consisting of 36 2kx4k, thinned, backside-illuminated 
CCDs that are each read out by two amplifiers. The mosaic has an unbinned pixel scale 
of O'.'OS which allows for a well sampled point-spread-function (PSF) even under the best 
seeing conditions. To decrease the read-out time we used 2x2 binning with the gain set so 
that the pixel sensitivity became non-linear before the analog-to-digital conversion threshold 
of 65,536 counts. Because of the fine sampling and the relatively deep pixel wells, one can 
collect 2 X 10^ photons in 1" seeing from a single star prior to saturation, setting the photon 
limit on the precision in a single exposure to ~ 0.25 mmag. 

The primary time-series photometric observations consist of ~ 4000 high quality images 
obtained over twenty four nights (including eight half nights) between December 21, 2005 and 
January 21, 2006. We obtained light curves for stars with 14.5 ^ r < 23 using a reduction 
pipelin e b ased on the image subtraction technique and software due to lAlard fc Lupton 
(119981 ) and lAlardI (120001 ). The resulting light curves were passed through the processing and 
transit detection pipeline that we describe in the following section. W e used the Daophot 



2 and Allstar PSF fitting programs and the Daogrow program (IStetson I Il987l . Il990 
19921 ) to obtain the g, r, and i single-epoch photometry. 



As desc ribed in Pap e r I w e also take BV photometry for stars in the field of this 
cluster from iKalirai et al.l (120011 ). Ks photometry from 2MASS (jSkrutskie et al.ll2006l ) and 
we tr ansform our ri photome try to Ic using a transformation based on the Ic photometry 
from iNilakshi fc Saeail J^OoJ • 



In addition to the photometry, we also obtained high-resolu tion spectroscopy of 127 st ars 
using the Hectochelle multi-fiber, high- dispersion spectrograph (jSzentgyorgyi et al. 1119981 ) on 
the MMT. The spectra were obtained on four separate nights between February 23, 2007 
and March 12, 2007 and were used to measure T^ff, [Fe/H], vsini and the radial velocity 
(RV) via cross-co rrelation agai nst a grid of model stellar spectra computed using ATLAS 
9 and SYNTHE jKurucj 119931 ). The classification procedure was developed by Meibom et 
al. (2008, in preparation), and made use of the xcsao routine in the iRAllfl rvsao package 



•^Iraf is distributed by the National Optical Astronomy Observatories, which is operated by the Associa- 



- 7- 



( iKurtz &: Minklll998l ) to perform cross-correlation. We use these spectra to provide stellar 
parameters and radial velocity measurements for several of the host stars to the candidate 
transiting planets. 



3. Transit Selection Pipeline 



In this section we describe the pipelines used to remove systematic variations from 
the light curves and select candidate transiting planets. We use these pipelines first to 
identify candidates (§ HI) and then to calibrate the detection efficiency by conducting transit 
injection and recovery simulations (§E])- The pipeline includes several steps: post-processing 
the light curves, devising transit selection criteria, and applying the selection criteria to the 
data to find candidates. We have devised three distinct pipelines, which we refer to as 
selection criteria sets 1-3, that differ in the number of post-processing routines applied and 
the manner in which candidates are selected. We first discuss the post-processing routines 
and the precision of the resulting post-processed light curves, we then discuss the selection 
criteria. 



3.1. Light Curve Post-Processing 

The raw light curves returned from the image subtraction procedure exhibit substantial 
scatter due to instrumental artifacts as well as astrophysical variations. Before searching 
these light curves for low-amplitude transit signals we must take steps to reduce the time- 
correlated noise. This process may reduce the sensitivity to high S /N planets, so for selection 
criteria set 3 we only apply a limited post-processing routine. Our post-processing pipeline 
consists of the following steps: 

1. (selection criteria sets 1-3) We first clip points from each light curve that are more 
than 5 standard deviations from the mean magnitude. We perform two iterations of 
this procedure. 

2. (selection criteria sets 1-3) For each image i, we determine /j, the fraction of light 
curves for which image i is at least a three sigma outlier. Images with /j > /c, 
the cutoff fraction, are removed from all the light curves. This process is performed 
independently for each of the 36 chips. We choose fc for each chip based on a visual 



tion of Universities for Research in Astronomy, Inc., under agreement with the National Science Foundation. 
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inspection of the histogram of fi values, we use values that range from 3% to 5%. This 
process typically removes ~ 500 images from the light curves. 

3. (selection criteria sets 1 and 2 only) We then remove 1 sidereal day, or 0.9972696 
day period signals from the light curves. This is done to remove artifacts due to, 
for example, rotating diffraction spikes that have a period of exactly 1 sidereal day. 
Figure [1] shows an example of a light curve that exhibits brightenings at a period of 
1 sidereal day as diffraction spikes from nearby bright stars sweep over the star. The 
diffraction spikes rotate in the image due to the need to rotate the camera with respect 
to the secondary mirror supports to keep stars on the same pixels throughout the 
night on the Alt-Az MMT. The number, width, depth, and shape of the brightenings 
seen in the light curves vary from star to star, moreover additional artifacts such as 
color dependent atmospheric extinction will also give rise to apparent variability with a 
period of 1 sidereal day. We remove these signals by binning the light curves in phase, 
using 200 bins, and then adjusting the points in a bin by an offset so that the average 
of the bin is equal to the average of the light curve. 

4. (selection criteria sets 1 and 2 only) As discussed in Paper III, we found that ~ 1/3 
of the probable cluster members that we observed show quasi-periodic variations with 
amplitudes of ~ 1% and periods ranging from 0.3 to 15 days. These variations are 
due to the presence of large spots on the surfaces of these relatively young (550 Myr), 
rapidly rotating stars. While it may be possible to identify a ~ 1% transit on top of 
such a signal, it would be increasingly difficult to identify shallower transits w ithout 



taking steps to remove these variations. Using the Lomb-Scargle algorithm (ILomb 



19761 : IScarglelll982l : iPress fc Rybickilll989l ) we identify the period P of the best-fit sine 



curve to each light curve. We then fit to the post-step-3 light curve a signal of the form 

m = Ao + Ai sin(2Pt/27r + 0i) + A2 sin(Pt/27r + ^2) + A3 sin(0.5Pt/27r + ^3) (1) 

where the Ai and 0j are free parameters, and we calculate Ax^^j.^^, the reduction in 
after subtracting the best-fit model from the light curve. Note that we adopt the 
convention that a more negative value of Ax^ indicates a better fit of the model to 
the data. This process increases the sensitivity to shallow transits but decreases the 
sensitivity to deep transits. We therefore fit a box-car transit signal to the post-step-3 
light curve , without subtracting the sine c urve, phased at period P using the BLS 



algorithm (IKovacs. Zucker and Mazehll2002l ) and calculate AxIlshx' reduction in 



after subtracting the box-car transit model from the light curve. We subtract equa- 
tion [1] from light curves with Axjii^j-m < Axblshx' i-^- "^^ '^'^^y ^PPly the correction to 
light curves for which a harmonic series fits the phased light curve better than a transit 
signal fits it. To test this technique we simulate light curves including variability due to 
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spots as well as transits. We use the iDorrenI (119871 ) spot model and the lMandel fc Agol 



(I2OO2I ) transit model. Figure [2] compares the results of a full BLS search on light curves 
with deep transits relative to spots and vice versa for the two cases of removing and 
not removing a harmonic signal before running BLS. We see that removing a harmonic 
series from light curves with Axiii^j-m < ^Xbls fix yields BLS results that are consistent 
with the injected transit signal, while not removing the harmonic signal from light 
curves with Axnarm > ^Xbls fix yields better results than removing the signal. 

5. (selection criteria sets 1 and 2 only) Finally we attempt to remove any remaining 
instrumental or weather related trends from the post-step-4 l ight curves using the 
Trend Filtering Algorithm (TFA iKovacs. Bakos fc NovesI boosh . This algorithm lin- 



early decorrelates each light curve against a representative sample of other light curves. 
The trend list for each light curve consists of the other stars on the chip with more 
than 2500 points, root-mean-square (RMS) < 0.1 mag and that are well outside the 
photometric radius of the star in question. There are typically ~ 250 stars in the trend 
list for each chip. 



3.1.1. Light Curve Precision 

In figure [3] we show the RMS as a function of magnitude for the stars after stages 2 and 
5 in the light curve processing pipeline. We plot both the point-to-point RMS as well as the 
RMS after binning the light curves by 2 hours in time. 



As discussed by lPont. Zucker and Queloz I (120061 ). the limiting factor for detecting low- 



amplitude transits around bright stars for many transit surveys is the presence of time- 
correlated systematic variations in the light curves (called red noise, noise that is uncorrelated 
in time will be referred to as white noise). We can estimate the degree of red noise in our 
light curves by fitting a noise model to the median RMS-r relation. We use a model of the 
form: 

where / = iQ-'^-^i™--^) is the effective number of photo-electrons per image for a source of 
magnitude m with zero-point z, s = iQ-o-^i"^s-z) jg ^.j^^ effective number of photo-electrons 
per image due to the sky that contaminate a given source and corresponds to a magnitude 
rris, and cr^ is the effective red noise in magnitudes. Note that z depends on the aperture 
and efficiencies of the optics and detector, the exposure time, the atmospheric extinction, 
and the photometric procedure. Also note that is the magnitude of a source which has 
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Q.^ G-S 

Fig. 1. — Example of a star with a light curve that shows repeated brightenings due to 
diffraction spikes from a neighboring bright star. The upper two panels show the star in 
question on a 0^5 x 0^5 FOV, when a diffraction spike is over the star (left) and when no 
spike is over the star (right). The two images were taken 100 minutes apart. The bottom 
panel shows the light curve of the star phased at 1 sidereal day, note the two periodic 
brightenings caused by the rotation of diffraction spikes in the images over the course of 
each night. 
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Fig. 2. — Transit recovery results for model light curves including both transits and photo- 
metric variations due to starspots. (Top) Comparison between the ratio of the model transit 
strength to the model spot strength, measured as a contribution to for each simulated 
signal, and the ratio of the best boxcar fit to the best harmonic series fit at the peak L — S 
period. We divide the data into two classes, those with ^Xb-lsmI ^^^avm ^ ^ (^•®" ^ transit 
model fits the light curve better than a spot model), and those with ^')(m.s,?Lyil ^'xkax-ai < ^ 
(i.e. a spot model fits the fight curve better than a transit model). (Bottom) Comparison 
between the best-fit BLS transit model, using a full period search, and the strength of the 
injected transit. We show the results separately for simulations with Ax^lq g^/Ax^arm > 1 
(left) and simulations with Ax|ls fix/AXnarm < 1 (right). In each plot we show the results for 
both removing and not removing a best-fit harmonic series from the light curve before run- 
ning the full BLS search. The solid lines show the ideal relation where Ax|ls fuU — AXtransit- 
When the transit model fits the fight curve better than a harmonic series (left) the harmonic 
series should not be removed from the light curve, when the harmonic series fits the light 
curve better (right) it should be removed before searching for transits. 
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Point -to- Point Pr^eiBion 2fcr-Bt^;iPi«iT^ialoii 




Fig. 3. — Light curve RMS vs. r magnitude for ~ 22, 000 stars after completing steps 1-2 
(top panels) and steps 1-5 (bottom panels) in the light curve processing pipeline. We show 
both the point-to-point RMS (left) and the RMS after binning the light curves by 2 hours 
in time (right), which is approximately the time-scale of a transit by an HJ. The solid lines 
show the median relation for the data plotted in each panel, while the dashed lines in the 
bottom two panels show the median relations from the top two panels. The dotted lines 
show the expected median relation for the binned hght curves if the RMS for each hght 
curve were reduced by N^^'^, where is the average number of points in a time- bin. 
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a flux equal to that of the total sky background through the effective photometric aperture. 
We bin the light curves on time-scales of 5, 10, 30, 60, 120, 180, 240, 1440, 2880, and 7200 
minutes and calculate the median RMS-r relation for each binning. We then fit equation [2] 
simultaneously to all ten relations using the free parameters zt-, and ar^x where and 
(7f.,T depend on the time-scale, while is independent of the time-scale. Table [1] gives the 
parameters from the model fit to the median RMS-r relations. We also list in table [1] the 
r-magnitude of a source for which the red noise on a given time-scale is equal to the Poisson 



2.5 V/t + St 



In(lO) /t 

As expected, ct^^t decreases with increasing T. At the time-scales relavent for a transit (1-3 
hours), the red noise is ~ 0.9 mmag. We also note that the values of zt are consistent 
with the expected values given the gain of the detector and the measured zero-point for the 
reference image, assuming the effective exposure time for the "average" image in the light 
curve is ~ 1.0 minute. 

An alternative method to determine the degree of red noise in the data is to examine 
the auto-correlation function of the light curves. We find that stars fainter than r ~ 18.0 are 
effectively uncorrelated in time while brighter stars appear to be uncorrelated on time-scales 
longer than ~ 200 minutes. 



3.2. Transit Selection 

To identify the best-fit transit signal in each light curve we use the BLS algorithm. We 
apply two slight modifications to the algorithm: 

1. Rather than searching over a fixed range of fractional transit width q = t/P values 
(where r is the duration of the transit and P is the orbital period), we allow the q 
range to vary with the trial frequency. We do this by assuming a stellar radius range of 
Rmin to Rmax and taking qmin,max = 0.076 {Rmin^maxff^ whcre / is the trial frequency 
in days~^, and assuming that M/Mq R/Rq for lower main sequence stars. 

2. While the post-processing described above substantially reduces the systematic vari- 
ations in the light curves, some variations remain (see fig. [3]). These variations give 
rise to increased power at low frequencies in the BLS periodogram of each light curve. 
As a result, long period transits should be treated with greater caution than short 
period transits. To account for this we subtract a mean filtered periodogram from the 
raw periodogram for each light curve before selecting the peak frequency. To mean 
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filter the periodogram we replace each point in the periodogram by the mean value 
of the 200 points closest in frequency to the given point after applying an iterative 
3cr clipping. Testing this modification on light curves with injected transits shows a 
very slight improvement (< 1%) in the fraction of light curves for which the recovered 
period agrees with the injected period. 

We conduct both a high resolution and a low resolution BLS search on the light curves. 
The low resolution search is needed for the transit recovery simulations (§ IS.l.ip to finish in 
a reasonable amount of time. For the high resolution search we examine 20, 000 frequency 
values over a period range from 0.2 to 5.0 days using 200 phase bins at each trial frequency, 
using Rmin = OARq and Rmax = 1.3i?0. For the low resolution search we examine 3072 fre- 
quency values over the same period range, using the same number of phase bins. In searching 
for candidates with the low resolution search we use identical parameters to those used in 
conducting the transit recovery simulations (§ IS.l.ip . We use slightly different parameters 
for possible cluster members and for field stars (see the selections of these stars in § 15.1.21 
and § 15. 2p . For candidate cluster members we take Rmin = O.li?© and Rmax = 1.1 x Rphot 
where Rphot is the radius of the star estimated from its photometry assuming that it is a 
cluster member. For the field stars we take Rmin = O.IRq and Rmax = 1.3i?0. 



As discussed by B06, a simple selection on the Sigr ial Residue (equation 5 inlKovacs. Zucker and Mazeh 



20021 ) or the Signal Detection Efficiency (equation 6 in lKovacs. Zucker and Mazelj|2002l ) gen- 
erally yields a large number of false positive detections that must then be removed by eye. 
However, because it is very difficult to accurately include by eye selections in a Monte Carlo 
simulation of the transit selection process, it is necessary to devise automatic selection cri- 
teria that minimize the number of false positive detections to calculate a robust limit on the 
frequency of stars with planets. We discuss the three distinct selection criteria that we have 
devised in turn. 



3.2.1. Selection Criteria Set 1 

To devise the first set of selection criteria we inject fake transit signals into the light 
curves of 1366 stars that lie near the cluster main sequence on a CMD and choose criteria 
that maximize the selection of high signal-to-noise simulated transits, while minimizing the 
selection of false positives. The simulated transits have radii of 0.35, 0.71 and 1.0 Rj, and 
inclinations of 90°. For each star/radius we simulate ten transits with periods ranging from 
0.5 to 5.0 days and random phases. We use the relations between r-magnitude and mass, 
and r-magnitude and radius for the cluster that were determined in Paper I to estimate the 
mass and radius of each star. The light curves including injected transits are passed through 
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the pipeline described in the previous subsection before running the BLS algorithm on them. 

To calibrate the selection criteria we first must define the set of injected transits which 
we consider to be recoverable, we then adopt selection criteria that maximize the selection of 
these light curves while minimizing the selection of all other light curves. Note that we use 
information from the injected transits to define the recoverable sample, while we cannot use 
any of this information when defining the selection criteria. We consider the transit to be 
recoverable if Axbls ~ ^Xblso < —100, where Axb^s reduction in for the best-fit 

BLS model to the light curve with the injected transit, and Axblso ^^S^^ curve 

prior to injecting the transit. When a light curve satisfies this criterion the best-fit BLS 
model is strongly influenced by the transit signal. This can happen even if the identified 
period does not match the injected period, or a harmonic of the injected period, however, in 
these cases, detailed follow-up may eventually yield the correct period. Note that for transit 
selection criteria sets 2 and 3 we adopt the more conservative recoverability criterion that 
the recovered period must match the injected period, or one of its harmonics, to within 10%. 

We find that the following selection criteria accurately distinguishes between the recov- 
erable and non-recoverable transits: 



1. Following BOG we use BLS to identify both the best-fit transit signal and the best-fit 
inverse transit signal for each light curve. For the best-fit tra nsit signal we calculate 



the signal-to-pink- noise ratio (jPont. Zucker and Queloz II2006I ) via: 



SN^ = (4) 

where 8 is the depth of the transit, nt is the number of points in the transit, Nt is the 
number of distinct transits sampled, cr^ is the white noise, and is the red noise at 
the time-scale of the transit. To calculate the white noise for a light curve we subtract 
the best-fit BLS model from the light curve and set cx^ equal to the standard deviation 
of the residual. To calculate the red noise we bin the residual light curve in time 
with a bin-size equal to the duration of the transit and set Or equal to its standard 
deviation. Note that while this technique provides a convenient method to determine 
an individual red and white noise estimate for each light curve, it overestimates the 
noise for an uncorrelated signal by a factor of ~ -\/2. We select candidates that have 
SN > 10.0. Figure m shows this selection. 

2. As shown in figure HI some of the transit signals that are not expected to be detectable 
pass the selection from the previous step. We can further reduce these potential false 
alarms with negligible loss of detectable transits by the following selection. Let Axbls 
refer to the reduction in for the best-fit transit signal and Ax'bls2 refer to the 
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reduction in for the second non-inverse transit peak in the BLS spectrum. The hght 
curves passing selection 1 are then selected if they have P > 2.0 days and Axbls " 
Ax%LS,2 < -60 or P < 2.0 days and AxIls - ^xIls,2 > -lO.gSP^ 'is shown in 
figure O 

When the above selections are applied to the actual data, more than 100 light curves 
are selected. Many of these are obvious false positives which can be eliminated by applying 
a few additional selections. 

3. We reject stars with fewer than 1000 points in their light curves, most of these are near 
saturation and show strong systematic variations even after applying TEA. 

4. Many of the false positives are faint stars located near much brighter stars, these can 
be rejected by only considering stars with an average instrumental magnitude brighter 
than 17.0 (corresponding roughly to r < 22.0). 

5. Other false positives come from very noisy light curves that can be rejected by requiring 
the standard deviation of the residual light curve after subtracting the best-fit box-car 
transit model be less than 0.1 mag. 

6. Many false positives show an anomalous faint set of points that occur on only one 
night. Following B06 we require that the fraction of Ax'bls ^^^^ comes from one night 
must be less than 0.8. 

7. We reject light curves for which the best-fit box-car transit has a period between 0.99 
and 1.02 days or less than 0.4 days. 

8. Finally we require that psr degree of freedom of the light curve residual after 
subtracting the best-fit box-car transit model must be less than 2.5 times the expected 

per degree of freedom. The expected per degree of freedom is taken to be the 
median value as a function of magnitude for the full ensemble of light curves (figure [H]). 
Note that the values of x^ per degree of freedom for stars between 18 < r < 23 approach 
~ 0.5 rather than ~ 1.0. This is due to a bug in our differential photometry routine 
whereby the formal differential flux uncertainties returned are not set to the flux scale of 
the reference image. We note that a similar bug is present in the differential photometry 
routine of the ISIS 2.1 package on which our photometry routine is based. The flux 
scale for the reference image that we use is set to that of a 30 second exposure taken in 
good seeing conditions. The formal errors are overestimated for longer exposures that 
were taken in poor seeing conditions. We identifled this bug after completing most 
of the analysis presented in this paper, we do not, however, expect it to signiflcantly 
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Fig. 4. — Histograms of the signal pink noise ratio of light curves with injected transit signals. 
We show separate histograms for injected transits that should be recoverable {^Xbls ~ 
^Xbls,o < —100) and for those that may be unrecoverable. Light curves with SN > 10.0 
are selected. 
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Fig. 5. — The difference in reduction from the best-fit box-car transit model and the 
second best-fit model vs. the period of the best-fit model for light curves that pass the 
selection in figure HI Light curves below the dashed line are rejected. 
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affect the statistical results presented here. Note that since the photometry that we 
are using is dominated by red noise, rather than white noise, the formal photometric 
errors, which assume Gaussian, uncorrelated noise, do not accurately describe the 
uncertainties in the data. We have therefore used a cut on the signal-to-pink noise, 
which is determined empirically from the light curve, rather than relying on the formal 
uncertainties. The formal uncertainties do affect the relative weightings of points, 
modifying them may thus yield slight changes in the TFA fitting procedure as well 
as in the periods that are identified by BLS. If the formally correct weighting scheme 
had been adopted, the sample of selected candidates may have been slightly different, 
however, this is no different from making minor changes to the rather arbitrary selection 
criteria. For our statistical conclusions regarding the fraction of stars bearing planets, 
what is important is that we apply the same selections and weighting scheme to the 
transit injection simulations as we do to the actual data and that the post-selection 
vetting procedure that we apply to the actual data would not reject any of the injected 
candidates. Making a uniform change to the weighting scheme will affect the selection 
of real candidates and injected transits in the same manner, so the transit detection 
efficiency measured by the injection simulations will accurately reflect the detection 
efficiency for the survey. 



3.2.2. Selection Criteria Set 2 



The above selection criteria were developed by manually inspecting the results of transit 
recovery simulations and defining cuts which appeared, by eye, to distinguish between suc- 
cessful recoveries and non-recoveries. While this method has the advantage that the criteria 
can be enumerated and easily visualized, it has the disadvantage that the cuts are fairly sub- 
jective and there is no guarantee that the selections optimally distinguish between recoveries 
and non-recoveries. 

To complement the above selection criteria, we have devised a set of selections using the 
Support Vector Machine (SVM) classification algorithm which has the ad vantage of bei ng 
much less subjective, but the disadvantage of being difficult to visualize (IVapnikI Il995l : a 
discussion of the algorithm can also be found in lPress et al.l 120071 which we summarize here). 



This algorithm takes as input a set of training data which consists of m points, (xi,?/j), 
where Xj is an n-dimensional vector of measureable parameters describing object i (called 
features), and Ui is either +1 or —1 and is used to divide the data into two classes. The 
algorithm then searches for a real- valued function /(x) such that /(xj) > for i/i > and 
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Fig. 6. — psr degree of freedom of the light curves without injected transits after sub- 
tracting the best-fit box-car transit model. The solid line shows the median per degree of 
freedom as a function of magnitude. Light curves above the dashed line (2.5 times the solid 
line) are rejected as transit candidates. See the text for a discussion of why ^ per degree of 
freedom for stars between 18 < r < 23 approaches ~ 0.5 rather than ~ 1.0. 
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/(xj) < for Hi < 0. The function is assumed to take the form: 

/(x) = W ■ 0(x) + B (5) 

where (f){'x.) is a fixed n-to-A^-dimensional transformation with typically larger than n, and 
W is an A^-dimensional vector. To optimize / one looks for a vector W that is normal to a 
hyperplane that separates the 0(xj) with i/i = +1 from the 0(xj) with i/i = —1 and that has 
a maximum perpendicular distance from the points nearest to it (called support vectors). In 
general, it is not always possible to find such a hyperplane, so instead one seeks to optimize 
/ by minimizing: 

^WW + Aj^S, (6) 

i 

subject to the constraint 

> 0, (7) 
Z/i/(x,;) > 1 - i= l,...,m 

where Sj are free parameters which allow for discrepancies between the model and actual 
classes, and A is a fixed regularization parameter used to control the tradeoff between accu- 
rately classifying the training data and maximizing the perpendicular distance between the 
support vectors and the hyperplane. In practice, the problem is recast in Lagrangian formal- 
ism so that one specifies a kernel matrix with the property Kij = K{yii,yij) = 0(xi) ■ 0(xj) 
rather than the transformation 0(x). 

To apply the SVM algorithm to ou r prob lem of devising transit selection criteria we use 
the SVMlight packag^ (Joachims 1999 . 2002 ). Before using the algorithm we apply a set of 



simple cuts which we found to be necessary to minimize the number of false positives: 



1. Reject stars with an average instrumental magnitude fainter than 17.0 (corresponding 
roughly to r > 22.0). 

2. Reject stars that have RMS > 0.1 mag, where the RMS here is the standard deviation 
of the residual light curve after subtracting the best-fit box-car transit model. 

3. Reject stars with SN < 10.0, where SN is the signal-to-pink noise given by eqn. |H 
Note for training the algorithm we use a less restrictive cut of SN < 9.0. 

4. Reject light curves for which the best-fit box-car transit has a period between 0.99 and 
1.02 days or less than 0.4 days. 



''C source code for SVMlight is freely available at 



http: / / svmlight.joacliims.org 
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We train the algorithm on the simulated transit data described in § I5.1.1[ We use only ~ 5000 
of the ~ 250 million simulated transits to allow the algorithm to converge in a reasonable 
amount of time. We take Ui = +1 for simulated transits that have 0.95 < Precover/ Pinject < 
1.05 where Precover and Pinject are the recovered and injected transit periods respectively, and 
l/i = —1 for all other simulated transits. There are 15 features in the Xj vectors including: the 
recovered period, the fractional transit width, the transit depth, SN, the white noise, the red 

noise, AxL^, ^xIls/^xIlsmv^ (^xIl5,2 - ^xIls)/^xIls^ the fraction of AxIls that 
comes from one night, the number of points in transit, the number of points observed less 
than T minutes prior to transit ingress (r is the duration of the transit), the number points 
observed less than r minutes after transit egress, the number of distinct transits observed, 
and the ratio of P^r degree of freedom of the light curve residual after subtracting the 
best-fit box-car transit model to the expected per degree of freedom. When applying the 
algorithm we consider a transit to be recovered if the estimated value of y is greater than 
0.1 as we found a number of false positives for which the algorithm returned values between 
0.0 and 0.1. 



3.2.3. Selection Criteria Set 3 

As described in § [HI when the first two sets of selection criteria are applied to transit 
recovery simulations the detection efficiency for 1.5i?j and larger planets drops unexpectedly 
for relatively bright stars. These large radius planets yield deep transit signals in the light 
curves of bright stars which are then distorted by steps 3-5 of the processing routines de- 
scribed in § 13. 1[ To recover these large radii planets we have devised a third set of selection 
criteria. For this set of criteria we run BLS on light curves processed through steps 1 and 2 
of § 13.11 (i.e. we do not remove one sidereal day period signals from the light curves, remove 
harmonic signals from the light curves, or apply TEA to the light curves). We then apply 
the SVM classification algorithm to the BLS results using the same training scheme and set 
of features as in § I3.2.2[ We train only on 1.0, 1.5 and 2.0Rj simulations. 

4. Transit Candidates 

Out of a total of 10,899 light curves for sources that were detected in g, r and i, have 
an average instrumental magnitude brighter than 17 (r < 22) and more than 1000 points, 
we select 16 transit candidates. Table [2] fists the candidates, their coordinates, photometry 
and BLS parameters and the selection criteria/resolutions which selected them. Phased 
light curves for the candidates are displayed in figure [71 As we discuss below, 12 of these 
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candidates can be rejected as eclipsing binary stars by their large transit depths or by noting 
that the primary and secondary transit depths are unequal. Two of the candidates detected 
on the unprocessed light curves result from artifacts in the data which are removed by the 
processing routines. One of the candidates is a blend with a nearby deep eclipsing binary. 
One candidate Jupiter-sized transiting planet remains. 

As seen in figure El several of the candidates (80014, 110021, 120050, 160017, and 
170100) lie close enough to the cluster main sequence on the g — r and g — i CMDs to be 
selected as potential cluster members in § I5.H however all of these candidates can be rejected 
as either eclipsing binaries or a blend with an eclipsing binary (see § 14.11) . We can therefore 
only place an upper limit on the planet occurrence frequency of cluster members. 



4.1. Discussion of Individual Transit Candidates 

30137 - This source is an eclipsing binary as evidenced by the secondary eclipse having 
an unequal depth to the primary eclipse when phased at twice the period shown in 
figure [71 This source is V1187 in the catalog of variables presented in Paper II. We 
reject this candidate. 

60161 - The transit-like variation in this light curve is completely removed by TFA so 
this source only passes selection criteria set 3. The transit feature is due primarily to 
2 nights. The images on these nights show a gradient in the background counts near 
the corner of Chip 6 where this star is located. The sense of the gradient is that points 
are fainter in the corner than near the center of the image. The use of a single master 
flat field image for the entire run appears to have failed for these two nights, the result 
is that many of the stars in the corner of Chip 6 show dips in their raw light curves on 
the two nights in question. This transit feature therefore appears to be spurious, and 
we reject it as such. We note that this is the only candidate that was selected by the 
high resolution search but not by the low resolution search. 

70127 - This source is a promising transit candidate around a field star. Cross- 
correlation analysis of the Hectochelle spectrum for the source yields Te// ~ 5000 K 
though the uncertainty is very high due to the faintness of the source. Note that 
there is a nearby source to 70127 that may also contaminate this spectrum. Combin- 
ing the -B, V and Ic photometry for the source, we estimate that [B — V)q ~ 1.0, 
and E[B — l^) ~ 0.16 so that the star has a radius of -R^ ~ O.75i?0 and a mass of 



Mi, ~ O.8M0. The light curve has iTingley fc SackettI (120051 ) parameters of rip = 0.7 



and r]* = 0.7, so it appears to be consistent with the transiting planet hypothesis. The 
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Phase 



Fig. 7. — Phased light curves for 16 candidate transiting planets selected by the pipeline 
discussed in § 13.21 See § 14.11 for a discussion of each candidate. The ID is taken from the 
photometric catalog presented in Paper I. We list for each candidate the period in days, 
the fractional transit width q and the transit depth in magnitudes returned by BLS. We 
note whether the displayed light curve has been processed through the full pipeline (TFA), 
through steps 1 and 2 only of § 13.11 (no detrend), or if no processing has been applied (raw). 
We show raw light curves for a handful of candidates where a-clipping removes some of the 
in-transit points. We also display a 9-bit flag to indicate which selection criteria sets and 
resolutions selected the candidate. From the left, the first three bits indicate if the candidate 
was selected by selection criteria sets 1, 2 and 3 respectively, using a high resolution BLS 
search. The middle three bits are for the low resolution BLS search applied to candidate 
cluster members, while the last three bits are for the low resolution BLS search applied to 
field stars. 
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Fig. 7. — continued. 
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B - V [Mag] V - [Mag] 

Fig. 8. — The location of 16 candidate transiting planets on g — r,g — i, B — V and V — Ic 
CMDs. The light points show all stars, the dark points show the transit candidates. Five of 
the candidates lie near enough to the cluster main sequence on the g — r and g — i CMDs to 
be selected as potential cluster members (see figure [9]). 



-27- 



source is located less than 2" from the edge of the chip, this appears to add some scatter 
to the light curve that is correlated with the seeing and removed by TFA. Also note 
that there is a 1.625 mag fainter source, 70181, that is located only O'.'S from 70127. 
The fainter source also shows a dip in its light curve, however the dip in flux appears 
to be greater for 70127 than 70181, and the centroid of 70181 on the residual image 
appears to shift slightly in phase with the transit, whereas 70127 does not. Both of 
these factors suggest that 70127 is the real variable. The 'shoulder' just before transit 
is an artifact of post-processing the fairl y high S/N light curve and does not appear 



in the unprocessed light curve. Fitting a iMandel fc Agoll (120021 ) transit model to the 



light curve w i th qu adratic limb darkening coefficients fixed to ui = 0.54 and U2 = 0.2 



from I Claret I (120041 ). which are appropriate for a ~ 5000 K dwarf star in the r filter, 
and fixing a/R^^ = 4.38 assuming the stellar mass, radius and the orbital period of 
0.77353 days, we find Rp ~ 1.0-Rj and sinz ~ 0.99. The expected RV amplitude for 
the star would be ii' ~ 220 m/s (Mp/Mj), where Mp is the mass of the planet. To 
rule out astrophysical false positives, such as an M-dwarf secondary, would require 
additional follow-up spectroscopy which would be challeng i ng giv en the faintness of 



the source (r ~ 19.6, V = 19.85). Note that ISahu et al.l (120061 ) achieved a formal 
RV precision of ~ 200 m/s for their candidate SWEEPS-11, which has a comparable 
magnitude (V = 19.83), so conducting spectroscopic follow-up for 70127 is not beyond 
the realm of possibility. 

80009 - This source appears to be an eclipsing binary, though there is not enough data 
to determine the period. There is a significant amount of scatter in the light curve that 
results from the source being nearly saturated. From the B, V, and Ic photometry we 
estimate the reddening to the source is E{B — V) ~ 0.18 so that the source would have 
{B — V)q ~ 0.65 and M ~ I.OMq. The primary eclipse depth appears to be deeper 
than the value returned by BLS (closer to ~ 0.1 mag), the secondary source would 
have to be significantly larger than 2.0Rj to create such a deep eclipse for this star. 
We reject this candidate. 

8OOI4 - This shallow transit-like signal is the result of a blend with the r ~ 18.3 mag 
eclipsing binary, V29, that is located ~ 3" away. We reject this candidate. 

90279 - This is an eclipsing binary with a period of 2.283 days. Phasing the light 
curve at twice the period shown in figure [7] reveals that the secondary and primary 
eclipses are slightly unequal in depth. It is VI 182 in the catalog of variables presented 
in Paper II. We reject this candidate. 



110021 - This source is a candidate cluster member with an estimated mass of ~ 
1.2M0. The secondary object must be significantly larger than 2.0Rj to yield an eclipse 
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deeper than 0.1 mag, the source is thus an echpsing binary. The source is V827 in the 
catalog of variables presented in Paper II. The period displayed in figure [7] and listed 
in table [2] is the period returned by BLS. When phased at a period of 6.7667 days the 
presence of a secondary eclipse becomes apparent, it is also apparent that the system 
has significant eccentricity, as well as an out of eclipse variation that phases with the 
proper orbital period. We reject this candidate. 

120050 - This is an eclipsing binary, note the dual eclipse depths at phase 0.0. When 
phased at twice the period shown in figure [7] it is clear that the light curve shows 
primary and secondary eclipses of differing depths. This source is V140 in the catalog 
of variables presented in Paper II. We reject this candidate. 

160017 - This source is an eclipsing binary. Its position on the CMD makes it a 
candidate cluster member with a mass of ~ 1.2Mq. For such a primary the deep 
eclipse could not be caused by a planetary-sized companion. There is an out of eclipse 
variation that phases at 4.0989 days, however the eclipses do not phase at this period. 
There is not enough data to accurately determine the orbital period. In figure [7] we 
show the raw light curve phased at the period returned by BLS for the light curve 
processed through TFA. The parameters listed in table [2] for this system are for the 
raw light curve, the S/N value is well below 10.0 in this case. The cr-clipping routine 
clips much of the eclipse yielding a shallow transit that phases roughly at the period 
shown and yields a S/N value above 10.0. This source is V791 in the catalog of variables 
presented in Paper II. We reject this candidate. 

160311 - This is an eclipsing binary, there is a shallow secondary eclipse with a depth 
of ~ 0.05 mag visible at phase 0.5. This source is V716 in the catalog of variables 
presented in Paper II. We reject this candidate. 

170049 - The transit feature seen in the noisy unprocessed light curve is completely 
eliminated by TFA. The source lies within 2" of the edge of the frame and the variations 
in the light curve are strongly correlated with the image seeing. The transit feature is 
removed when the light curve is de-correlated against seeing. We reject this candidate 
as the variations appear to be spurious. 

170100 - This source is an eclipsing binary. When phased at twice the period displayed 
in figure [7] it is clear that the primary and secondary echpses are of unequal depth. 
The out of eclipse variations also phase at this period. This source is V1028 in the 
catalog of variables presented in Paper II. We reject this candidate. 

230082 - This source, V485, is most likely an F-M eclipsing binary. Note that when 
phased at the period recovered by BLS (as shown in figure [7]) there appears to be a 
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secondary eclipse at phase —0.5, when phased at a period of 1.6676 days, however, 
the putative secondary matches to the primary echpse. Also note that the clipping 
procedure removes some points from the bottom of the eclipse, the eclipse in the raw 
light curve is ~ 0.01 mag deeper. From the B, V and Iq photometry we estimate 
E{B - 1/) ~ 0.38, so {B - V)o ~ 0.41 which corresponds to a ~ l.SM©, and ~ IARq 
primary star. For a star of this size, a 0.034 mag transit requires that the companion 
have a radius greater than 2 Rj. 

270149 - This source, V457, is an eclipsing binary. Note the deep, and distinct, primary 
and secondary eclipses. We plot the raw light curve for this star in figure [71 the clipping 
routine removed many of the points in eclipse which caused BLS to identify the wrong 
period and underestimate the depth. 

280287 - This is a grazing eclipsing binary. When the raw light curve is phased at 
twice the period found by BLS it is clear that the primary and secondary minima are 
of unequal depths. This source is V226 in the catalog of variables presented in Paper 
II. We reject this candidate. 

330224 ^ This source, V141 in the catalog of variables, is an eclipsing binary given its 
depth. The raw light curve reveals a strong out of eclipse variation with a peak-to-peak 
amplitude of ~ 0.1 mag that phases at the orbital period. 

5. Transit Detection Efficiency Calculation 

To calculate our planet detection efficiency we follow the procedure described by B06; 
we summarize the procedure here and perform the calculation separately for cluster members 
and field stars. 

5.1. Cluster Members 

In figure M we show the selection of stars near the cluster main sequence on g — r and 
g — i CMDs. We select a total of 2475 stars that have an instrumental magnitude less than 
17, RMS < 0.1 mag and that have at least 1000 points in their light curves. Note that we 
expect ~ 1450 of these stars to be cluster members. We use these stars to place a limit on 
the planet occurrence frequency for the cluster. The total number of planets in the cluster 
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that we expect to detect with our survey is given by 



N 

Ndet = / Pdet,i 



8 



i=l 



where / is the planet occurrence frequency over a specified planet radius and orbital period 
range, Pdet,i is the probability of having detected such a planet for star i, and there are N 
candidate cluster members in the survey. For a binomial distribution, the detection of no 
planets is inconsistent at the ~ 95% level when N^et ^ 3. The 95% confidence upper limit 
on the planet occurrence frequency is then 



where P^,* is the probability of detecting a transit of a planet with orbital period P and 
radius Rp for star i if it has an orbital inclination that yields transits, Py.i is the probability 
that the planet has an inclination that yields transits, Pmem,i is the probability that star i is 
a cluster member, and d^p/ dRpdP is the joint probability distribution of Rp and P. 

There are four terms that contribute to equation [10] that must be determined separately. 
Assuming random orientations, the term P^,* is given analytically by Pt = (R^+Rp) / a, where 
R^ is the radius of the star and a is the semi-major axis for an orbital period P and stellar 
mass M^. For the term Pmem,i we take the photometric membership probability of the star, 
which we calculate by comparing the CMD of the cluster to a CMD of a field off the cluster. 
The term d^p/dRpdP must be treated as a prior, and the term P^^j is calculated via Monte 
Carlo simulation as discussed below. 



We follow the procedure described by B06 to calculate P^ i. This involves injecting 
limb-darkened transits into the light curves of potential cluster members and attempting to 
recover them. Transits are injected into the raw light curves, we then run each simulated 
light curve through the post-processing and transit selection routines described in § |3l Note 
that we only conduct simulations on light curves that are not selected as transit candidates. 




(9) 



The planet detection probability for star i is given by 




(10) 



5.1.1. Calculating P^^i 
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Fig. 9. — The selection of stars near the cluster main sequence, which are used in computing 
the planet detection probability, are plotted on g — r and g — i CMDs. The bottom two panels 
show the CMDs for the field on the cluster, while the top two panels show the CMDs for a 
field adjacent to the cluster with the same Galactic latitude (see Paper I). The off-cluster 
field is used in determining the membership probability as a function of magnitude (fig. [T0|) . 
Stars falling within the solid black lines on both CMDs are selected as stars near the cluster 
main sequence. 
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Since there are only a few candidates, not including their efficiences should not change our 
results significantly. 



To simulate transit light curves we use the iMandel fc Agoll (120021) analytic model with 
r-band quadratic limb-darkening coefficients from I Claret I (120041 ) . Note that we assume 
circular orbits for the short-period planets under consideration. We estimate the mass. 



radiu s and surface temperature of each star using the best-fit YREC isochrone (lAn et al. 



20071 . also see Paper I), we th en estirnate th e limb-darkening coefficients for each star via 
linear interpolation within the IClaret I (12004J ) grid, assuming [M / H] = 0.045 (Paper I) and 
Vturb = 2 km/ s. 

To determine the dependence of P^^i on the orbital period we inject the transits in period 
bins ranging from 0.4 to 5.0 days with a logarithmic step size of 0.022. For each period bin 
we inject Ntriai transits with periods distributed uniformly in logarithm over the bin, random 
phases and inclinations distributed uniformly in cos i over the range < cos i < (i?* -|- Rp) / a. 
Following BOG we estimate that the error in the recovery fraction / = Nrecover / Ntriai is given 
by 



trial • 



For each period bin we initially adopt Ntriai = 20 and continue simulating transits until 
aj < 0.05 using selection criteria set 1. We conduct simulations for planetary radii of 0.3, 
0.35, 0.4, 0.45, 0.5, 0.6, 0.7, 0.8, 0.9, 1.0, 1.5 and 2.0Rj. 



5.1.2. Calculating Pmem,i 

To calculate the membership probability for each star we use the luminosity functions 
of the cluster and Galactic field along a strip in the g — r and g — i CMDs enclosing the 
cluster main sequence (see figure [9], the determination of the luminosity functions is described 
in Paper I). This gives the membership probability as a function of r-magnitude only. In 
figure flOl we show the membership probability as a function of r. Since this method ignores 
color information in assigning a probability to stars, it tends to give too high a probability to 
stars lying away from the cluster main sequence, and too low a probability to stars lying close 
to the main sequence. If the transit detection probability depends only on the r magnitude 
of the star then this simplification should not bias the result. However, as shown in Paper 
II, cluster members have a higher probability of exhibiting photometric variability than field 
stars of the same magnitude, so we caution that it is possible that the transit detection 
probability may be lower on average for cluster members than for field stars. We consider 
this possibility in § 16. 3[ 
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Fig. 10. — Histogram of stars in the field of the cluster selected in figure [9] as a function of r 
magnitude (open histogram) compared to the histogram of cluster members (filled histogram) 
taken to be the difference between the histogram of the field on the cluster and the histogram 
of the field off the cluster. The solid line shows the membership probability (calculated as 
the ratio of the filled histogram to the open one) as a function of r. The histograms are taken 
from Paper I and include a correction for photometric incompleteness, but do not correct 
for spatial incompleteness. This figure is analogous to fig. 9 in B06. 
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5.2. Field Stars 

To calculate the detection efficiency for field stars we follow a procedure that is similar 
to what we use for cluster members. In this case, however, the mass and radius of a star 
cannot be determined simply from the magnitude of the star. Instead we use Galactic models 
to estimate the mass and radius for each star. 

First we select a sample of observed field stars that includes all stars with g, r, i, B and 
V photometry that have more than 1000 points in their light curves, have a light curve RMS 
that is less than 0.1 mag, an instrumental r-band magnitude less than 17 (corresponding 
roughly to r < 22) and were not selected as candidate cluster members. The sample of stars 
is shown on B — V and V — Ic CMDs in figure [TTl A total of 7824 stars are selected in this 
fashion. We then obtain simulated photometric observations of a 24' x 24' field centered a t 



the Galactic latitude/longitude of the cluster using the Besangon model (iRobin et al. II2003I ). 
We assume an interstellar extinction of 0.7 mag/kpc. We caution that this model is known 
to be unreliable a long certain lines of sight, in particular for low Galactic latitudes (e.g. see 



Ibata et al.l 120071 ). For each observed star in our sample we choose the simulated star that 



mmimizes 

{Vo - Vsf + {Bo - B,f + (/c,o - Ic,sf (12) 

where the o subscripts denote photometric measurements for observed stars, and the s sub- 
scripts denote photometric measurements for simulated stars. We then assign the mass and 
radius of the simulated star to the observed star. We reject 10 stars which match to white 
dwarfs since these all lie in an isolated region of the CMD. The majority of the stars in 
the sample (97%) match to dwarf stars (log (7 > 4.0, cgs) rather than sub-giants or giants. 
Figure [T^] shows the estimated radii of the field stars as a function of magnitude compared 
with the values for the cluster. Note that at fixed magnitude, the majority of field stars have 
larger radii then the cluster stars. As a result, we expect the transit probability to be larger, 
but the overall S/N to be smaller at fixed magnitude for the field stars. Thus for signals for 
which the S/N is much larger than the threshold, the detection efficiency is larger for field 
stars, but for planets near the S/N threshold, the detection efficiency will be smaller. Thus 
the minimum detectable planet radius will be smaller for cluster stars, all else being equal. 
For the less common foreground field dwarfs, the opposite is true. 

We determine the detection efficiency by conducting transit injection and recovery sim- 
ulations as described in § 15.1.11 We conduct the simulations for 1000 randomly selected field 
stars. For each star we calculate Pdet,i using eq. [TO] where we now take Pmem,i = 1- The 95% 
upper limit on the occurrence frequency, assuming no detections, can then be calculated us- 
ing eq. [9]with the sum in the dominator multiplied by 7814/1000 to scale from the simulated 
subset to the full sample. 
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Fig. 11. — B~V (left) and V — Ic (right) CMDs for stars selected as members of the Galactic 
field (top), for a simulated observation using the Besangon Galactic models (middle) and for 
a simulated observation using the Trilegal Galactic model (bottom). For plotting purposes 
we have added Gaussian noise of 0.05 mag to the V — Ic colors from the Besangon model. 
Stars not selected as cluster members that pass a number of cuts on magnitude, light curve 
RMS and the number of points in the light curve (see § 15.11) are selected as members of the 
Galactic field. The solid lines show approximately the location of stars selected as potential 
cluster members, these stars are selected on gri CMDs and are not included in the top two 
CMDs. To estimate the mass and radius for each star in the top two CMDs we take the 
values for the nearest star in the middle two CMDs. We use the Trilegal model to estimate 
the errors in the planet frequency upper limits that result from uncertainties in the Galactic 
model (§ 16.31) . Open stars in the top two plots show points that match to white dwarfs. We 
reiect these stars from our samnle. 
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Fig. 12. — The estimated stellar radius is plotted against the magnitude for Galactic field 
stars (solid points) when using the Besancon Galactic model (top) and the Trilegal model 
(bottom). The solid line shows the relation for the cluster. 
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6. Results 

In figure [13] we show several examples of transit injected light curves recovered in our 
simulations. Note that 1.5 and 1.0 Rj planets are easily recovered. Also note that the S/N 
for .4-R.7 planets is similar to that f or ~ l.ORj planets discovered by the field surveys (see, 
e.g. 



Pont. Zucker and Queloz II2006I ). 



In figure[T4]we plot the detection probability (eq. [TO]) as a function of magnitude for 0.5, 
1.0 and 1.5Rj simulations on candidate cluster members recovered using selection criteria 
set 2. We show the results for Extremely Hot Jupiter (EHJ, 0.4 < P < 1.0 days). Very Hot 
Jupiter (VHJ, 1.0 < P < 3.0 days) and Hot Jupiter (HJ, 3.0 < P < 5.0 days) period ranges. 
For the HJ period range we use an upper limit of 5.0 days rather than 9.0 days since we do 
not attempt to recover planets with periods longer than 5.0 days. For each period range and 
planet radius we integrate eq. [10] assuming a uniform probability distribution in log P and a 
Dirac-5 function distribution in R, i.e. 

72 

^ ocSiR-Ro)/P (13) 



dPdR 



where Pq is the planet radius under consideration. Note that the detection probability for 
1.5 Rj planets orbiting bright stars (15.0 ^ r < 20.0) drops relative to the probability 
for 1.0 Rj planets. There are two factors that contribute to this: steps 3-5 of the light 
curve processing routines described in § 13.11 distort the high signal-to-noise transits for these 
large planets often leading to out-of-transit variations in the processed light curves, and 
the discrepancy between the limb-darkened transit signal and the box-car model becomes 
significant for the high signal-to-noise transits. When the light curve processing routines 
are not used (selection set 3) the drop in the detection probability of large planets is less 
significant, however the sensitivity to smaller planets is reduced. 

In figure [15] we compare the detection probability for field stars to the probability for 
cluster members, we set Pmem = 1 for the cluster members in making this comparison. 
As expected, the field stars generally show slightly higher detection probability at fixed 
magnitude than the cluster members for signals with S /N well above the detection threshold 
due to the higher transit probability, but smaller for those signals near the threshold due to 
the larger radii (and so shallower transits). 

Figure [TB] shows the ~ 95% and 1 — a confidence upper limits on the planet occurrence 
frequency for cluster members and field stars as a function of orbital period using selection 
criteria set 2. We assume the orbital periods are uniformly distributed in logarithm within 
logarithmic period bins of size 0.022 and show the results for several different planet radii. 
In figure [T7] we show the ~ 95% confidence upper limits for cluster members as a function of 
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Fig. 13. — Example phased light curves of simulated transits that are successfully recovered 
using selection criteria set 2 (§ l3.2.2p . The simulated planets have radii of 1.5, 1.0, and 0.4 Rj 
while the stars have radii of 0.70, 0.47, and 0.83 Rq respectively. The injected periods are 
1.77, 1.44 and 1.12 days respectively. The line shows the best-fit box-car transit, we also 
list the signal-to-pink- noise ratio for each transit. Note that the examples shown correspond 
approximately to the median signal-to-pink-noise recoveries for each planet radius among all 
recovered simulations with periods between 1.0 and 3.0 days. The 1.5 and 1.0 Rj planets 
are easily recovered while the 0.4 Rj planet is marginally recovered. 
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Fig. 14. — The detection probability as a function of r-magnitude for transit injec- 
tion/recovery simulations of candidate cluster members. The simulations are recovered using 
selection criteria set 2. We show the results for the HJ, VHJ and EHJ period ranges for 

several planetary radii. The detection probability decreases toward fainter magnitude since 
the transit probability at fixed period decreases for smaller stars. Note that for very high 
signal-to-noise transits (1.5 Rj) the detection probability is lower than expected as a result 
of the trend-filtering routines. 
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Fig. 15. — The detection probability as a function of r-magnitudc for field stars (dark points) 
compared witli candidate cluster members (light points). We plot PcPt which corresponds 
to Pdet when Pmem ==1-0 for cluster members and is equal to P^et for field stars. 
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planet radius for the EHJ, VHJ and HJ period ranges. We compare the results for the three 
distinct selection criteria sets. The results are listed as well in table [31 We find that the SVM- 
based selection (selections 2 and 3) outperforms the non-SVM-based selection (selection 1). 
The upper limits set using SVM are as much as 1.7 times smaller than the upper limits set 
using the non-SVM selection. 

In figure [18] we show how the distribution of stellar masses to which we are sensitive to 
planets depends on planetary radius and period for the cluster candidates and field stars. 
Note that for 0.3-Rj planets the field star distribution is peaked toward smaller stars, whereas 
the cluster distribution is peaked toward larger stars. The small field stars to which we have 
sensitivity to 0.3-Rj planets are foreground stars, at the distance of the cluster these stars 
are too faint for us to detect Neptune-sized planets around them. For larger planets the field 
star sensitivity distribution is peaked around O.S-I.OMq, whereas the cluster distribution is 
broader and peaked at around O.SMq. 

6.1. Cluster Members 

For cluster members we place a 95% confidence upper limit on the frequency of stars 
with Jupiter-sized EHJ, VHJ and HJ planets of 1.1%, 2.7% and 8.3% respectively. Note that 
all of these limits, as well as the ones we discuss below, come from selection criteria set 2. 
For smaller planets the limits rise. For the EHJ period range we can place a limit of 25% on 
planets down to 0.35 Rj, which is roughly the size of Neptune. For the VHJ period range, 
we can place a limit of 44% on planets down to 0.45 Rj, while for the HJ period range we 
can place a limit of 49% on planets down to 0.6 Rj. 

6.2. Field Stars 

The sample of field stars provides more stringent upper limits on the planet occurrence 
frequency than the sample of cluster members. For this sample we can place a 95% confidence 
upper limit on the frequency of stars with Jupiter-sized EHJ, VHJ, and HJ planets of 0.3%, 
0.8% and 2.7% respectively. This assumes that the candidate 70127 is not a real planet, if 
70127 is real the frequency of ~ 1.0 fZj planets in the EHJ period range would be 0.002% < 
/ < 0.5% with 95% confidence. In principle this frequency is small enough that planets 
in this period range could have escaped detection in most other radial velocity and transit 
surveys. Some transit surveys have probed enough stars to be sensitive to planets with these 
frequencies, as some of them have looked at as many or more stars than we have, but these 
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Fig. 16. — The ~ 95% confidence upper limit (solid lines) on the fraction of cluster members 
(left) and field stars (right) with planets of a given radius as a function of orbital period. We 
show the results for 0.35, 0.5, 0.7 and 1.0 Rj planets. For field stars we show 0.3 Rj rather 
than 0.35 Rj. We also plot the 1 — a upper limits (dotted lines) for each of these planet 
radii. 
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Fig. 17. — The ~ 95% confidence upper limit on the fraction of cluster members (left) and 
field stars (right) with planets in a given period range as a function of planetary radius. We 
show the resuhs for the EHJ (0.4 < P < 1.0 days), VHJ (1.0 < P < 3.0 days) and HJ 
(3.0 < P < 5.0 days) ranges. We compare the results for each of the selection criteria sets 
discussed in § 13. 2[ 
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Fig. 18. — The Pdet weighted distribution of steUar masses to which we have sensitivity 
to detect planets. The curves have been normahzed to have unit integral over the range 
O.OQMq < M < 2.1Mq. The distributions are plotted for different planetary radii and 
period ranges, and are shown for cluster candidates (left) and field stars (right) separately 
Note that for 0.3i?j planets the field star distribution is peaked toward smaller stars, whereas 
the cluster distribution is peaked toward larger stars. The small field stars to which we have 
sensitivity to 0.3it!j planets are foreground stars, at the distance of the cluster these stars 
are too faint for us to detect Neptune-sized planets around them. For larger planets the field 
star sensitivity distribution is peaked at around O.S-I.OMq, whereas the cluster distribution 
is broader, but also peaked at around O.SMq. 
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surveys have gene r ally not been very sens itive to hosts with masses as small as this host (see 



Gould et al.l l2006l : iBeatty &: Gaudil l2008l ) . For the EH J period range we can place a limit 
of 16% on planets down to 0.3 for the VHJ range a limit of 8.8% on planets down to 
0.5 Rj and for the HJ range a limit of 47% on planets down to 0.5 Rj. Extrapolating the 
curves in figure [T71 we note that for the EH J period range we do have some sensitivity even 
to planets as small as ~ 2.5/?®. 



6.3. Errors on the Upper Limits 

There are a number of factors which may contribute errors to the upper limits. This 
includes uncertainties in from using a finite number of simulations, errors in Pmem from 
using a finite sample, and systematic errors from blends and binaries, from errors in the 
assumed Galactic model and from the possibility that cluster members are more likely to 
be variable than non-cluster members. BOG give a detailed discussion of how most of these 
errors can be estimated. 

The fractional uncertainties in the 95% upper limits for cluster members (o"/<g5%//<95%) 
due to using a finite number of simulations to determine range from ~ 1% for 0.3i?j 
simulations to ~ 0.04% for 1.0i?j EHJ simulations. For field stars the fractional uncertainties 
are all less than 1%. These errors are negligible compared to other sources of uncertainty. 

Following B06 we note that the fractional uncertainty in the 95% upper limit for candi- 
date cluster members due to uncertainties in Pmem is equal to the fractional uncertainty in 
the effective number of cluster members: a^^ ^^^/N^,^eff, where N^^^^ff = < Pmem >■ We 
find aN,^^^^.^./N,^,ff = 2%. 

We expect binarity to be a more significant effect than chance alignments. As we noted 
in Paper II, we expect ~ 1 chance alignment within O'.'l of two point sources brighter than 
our detection threshold in our entire field. The effect of binarity on a transit survey is not 
straightforward. While the blending of light from two stars will dilute the transit signal, the 
primary star will be slightly smaller than what is assumed when injecting transits. Moreover, 
if the precision of the light curve is good enough that one could still detect the transit for 
a given planetary radius and primary star radius if the transit were more than a factor 
of 2 shallower, then such a planet would be detectable orbiting the primary star for any 
mass ratio and it may also be detectable orbiting the secondary star. In this case the total 
number of stars to which one is sensitive to planets is grea ter than estimated when binarity 



is neglected (see iGould et al.ll2006l : iBeatty fc Gaudil l2008l . for further discussion about the 



effects of binarity). Binarity will only affect the detectability of a transit if the mass ratio is 
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high, B06 estimate that the requirement is g > 0.6, which they argue is true for only ~ 11% of 
dwarf stars when the binary fraction ~ 50%. The re is an indication th at the binary fraction 



along the main sequence in M37 is closer to 20% (iKalirai &: Tosill2004j ). which would reduce 
the number of high mass ratio binaries to ~ 4% assuming that the mass ratio distribution 
in the cluster is the same as for the Galactic field. 

To estimate how uncertainties in the inferred masses and radii of field stars that results 
from uncertainties in the Galactic model affect the planetary fre quency upper limit s, we 



recompute the upper limits using the Trilegal 1.2 Galactic model (IGirardi et al.ll2005l ). We 
match the sample of field stars to a simulated set of photometric observations generated 
with the Trilegal model following the same procedure that we used in § 15.21 to match to 
the Besangon model. The Trilegal model is shown in figures [TT] and [T2] and is generated 
assuming a local extinction law of 0.7 mag/kpc. With this model we find that 96.5% of stars 
in our field sample have logg > 4.0, which is comparable to the dwarf fraction found with 
the Besangon model. In § 15.21 we computed the transit detection probability for a sample 
of 1000 field stars by conducting transit injection/recovery simulations with stellar masses 
and radii adopted from the match to the Besangon model. Here we avoid the expensive 
task of conducting additional transit injection/recovery simulations by matching stars in the 
Trilegal-matched sample to appropriate stars in the Besangon-matched sample. We match 
each star, i, in the Trilegal-matched sample to the star, j, in the Besangon-matched sample 
that minimizes 

(0.21n(10)(r, - r,))^ + (2(i?, - Rj)/R,y (14) 

where r is the r-magnitude of the star, and R is its radius. We then set the transit detection 
probability for star i equal to the probability for star j determined from the simulations. 
This choice of weighting between magnitude and radius minimizes differences between the 
expected transit signal to noise ratio of the two stars. The upper limits on the planet 
occurence frequency are then computed with eq. [91 The results using selection criteria 2 
are listed in table |H We find that the fractional difference in the planet frequency upper 
limits from the two Galactic models is < 10%. The Trilegal model yields a systematically 
lower upper limit for the O.SRj and 0.5Rj planets but a systematically higher upper limit 
for larger planets. In a similar manner we recalculate the upper limits using the Besangon 
model generated with less extinction (0.5 mag/kpc). In this case we find that the upper 
limits for planets smaller than l.ORj are systematically smaller while the upper limits for 
larger planets are systematically larger. The fractional differences in the upper limits for 
0.7Rj and smaller planets range from ~ 5% to ~ 30%, while for larger planets they range 
from ~ 2% to ~ 20%. We adopt a fractional uncertainty of ~ 10% on the upper limits due 
to uncertainties in the Galactic model. 



To estimate the effect of variability on the upper limits for cluster members we compute 
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the upper limits under the extreme case where all candidate cluster members that are variable 
have Pjnem = 1 and other stars have Pmem = {N{r)Pmem,o{r) - N^ar{r))/{N{r) - N^ar{r)) 
where Pmem,o is the value of Pmem when variable and non-variable stars are weighted equally, 
N[r) is the number of candidate cluster stars in magnitude bin r and N^ar is the number 
of variable candidate cluster stars in magnitude bin r. The resulting 95% upper limits for 
selection criteria set 2 are given in table [51 For small radius planets the effect is significant, 
so that the upper limit on 0.35Rj EHJ planets, for instance, increases to 38% from 25%. For 
larger planets the effect of variability is less important, above l.ORj the fractional increase 
in the upper limit is less than 10%. 

Assuming that binarity and variability only increase the upper limits, the fractional 
uncertainty on the upper limits for 1.0-Rj planets orbiting candidate cluster members is 
~ +11%, —2%, for 0.35-Rj planets it is ~ +50%, —2%. For field stars we can neglect the 
uncertainty due to variability and the uncertainty on the membership probability, but we 
must include the uncertainty on the Galactic model, so that the fractional uncertainty on 
the upper limit for all planetary radii is ~ +15%, —10%. 



6.4. Comparison with Previous Results 



Our results can be compared both with previous results for transit surveys of open 
clusters, and with results from transit and RV surveys of the Galactic field. 

Our limits for the cluster on the frequency of EHJ, VHJ and HJ planets with R = 1.5Rj 
of 1.2%, 2.3% and 6.9% respectively are substantially better than the corresponding limits 
for NGC 1245 by B06 of 1.5%, 6.4% and 52%. The primary differences between the two 
surveys are that we observed approximately twice as many cluster members as B06 and we 
obtained significantly higher precision photometry at fixed stellar radius than B06. Note, 
however, that B06 use a period range of 3 — 9 days for the H J range whereas we use a range of 
3 — 5 days. Since, by selection, we have zero sensitivity to planets with 5 days < P < 9 days, 
extending our HJ range to the B06 range results in a limit of ~ 15% on the fraction of stars 
with planets o f radius 1.5 R.j (a ssuming an underlying distribution in P that is uniform 
in logarithm). iMiller et al.l (120081 ) have also conducted detailed Monte Carlo simulations of 
their open cluster transit survey and placed 95% upper limits of 14% and 45% on VHJ and 
HJ with radii of 1.5-Rj. 

We can estimate the frequency of planet s detected by RV surv eys in the EHJ, VHJ 
and HJ period ranges using the results from ICumming et al.l (l2008l . see figure 5 in their 
paper). Out of their sample of 585 FGKM stars, there are 7 planets with M > O.lMj 
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and 3 days < P < 5 days, 1 planet with M > O.lMj and 1 day < P < 3 days, and no 
planets detected with P < 1 day. This yields 95% confidence intervals for the occurrence 
frequencies of 0.48% < fnj < 2.5%, 0.0043% < fvHj < 0.95%, and /ehj < 0.51%. Our 
results are not directly comparable to the RV survey results for two reasons: first, our 
limits are based on planetary radius while the RV survey limits are based on planetary mass, 
second, the distribution of stellar masses and metallicities is not the same for the two surveys. 
Nonetheless, if we assume all extra-solar giant planets have radii of ~ l.OPj and that the 
planet occurrence frequency does not depend on stellar properties, then our upper limits of 
fsHj < 0.3%, fvHj < 0.8% and fnj < 2.7% for field stars are consistent with the measured 
RV frequencies. 



Gould et al.l (120061 ) used the results of the OGLE transit survey to determine the fre- 
quency of VHJ and HJ planets. Assuming planets follow a uniform distribution in radius 
between l.OPj and 1.25Pj, they found fvHj = 0.141°,^^% and fnj = 0.3llJ^:?^% with 90% 
confidence error bars. Furthermore, they placed a 95% confidence upper limit of < 1.7% on 
planets with perio ds between 1 and 5 days and radii distributed uniformly between 0.78 and 



0.97Pj. Similarly iFressineLaLl (120071) found fvHj = 0.18% and fnj = 0.29% for the OGLE 
survey, and they found that the fraction of stars bearing planets with periods less than 2 
days is 0.079]'^o;q4o% with 90% confidence error bars. Note that the frequency of HJ is lower 
for the OGLE survey than the RV surveys because the RV surveys tend to be biased towards 
higher metallicity stars whereas the OGLE survey is not. Our upper limits on the VHJ and 
HJ occurence frequencies for field stars are consistent with the frequencies determined by 
both of these groups. 



W e can also compare our results to the results from the SWEEPS survey. ISahu et al. 
(120061 ) find that ~ 0.4;1^q;2% of stars larger than 0.44 Mq host a Jupiter-sized planet with 



P < 4.2 days, assuming all 16 of their candidates are real, which is consist ent with the limits 



that we set. Focusing specifically on the EHJ period range, we note that ISahu et al.l (120061 ) 
found 5 candidate planets in this range and estimate that at least ~ 2 of them are likely 
to be real. Assuming all 16 candidates are real planets, and that 5 of them are EHJ, their 
results suggest that ~ 0.13% of stars host an EHJ. If, on the other hand, we estimate that 
only half of the candidates are real, and that 2 of them are EHJ, the frequency of EHJ would 
be ~ 0.05%. Note that we have assumed that the detection probability is constant over the 
entire period range. Accounting for the enhanced probability of detecting EHJ over VHJ 
and HJ may lower the frequency estimates of these planets by as much as a factor of ~ 5. 
The SWEEPS frequency of EHJ is below the 95% upper limit that we set assuming 70127 is 
not a real candidate. If 70127 is real, our EHJ frequency of 0.09]'^q;q77% (with la errorbars) 
would be consistent with the SWEEPS results. 
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6.5. Expected Planet Yield 

Based on the planet dist ribution determi n ed by other surveys we can estimate the 



expected yield of our survey. ICumming et al.l (120081 ) find that the distribution of planet 



masses down to ~ O.lMj is given by 



dlnMdlnP 



CWpf^ (15) 



where a = —0.31 ± 0.2 and f3 = 0.26 ± 0.1, M is in Jupiter masses, P is in days, and C 
is a constant. For notational simplicity, in the following discussion we do not write factors 
of Mj and Rj. We will adopt the mass dependence from this relation, but not the period 
dependence, since this relation does not account for the pile-up of planets at short periods. 
For the period dependence we adopt a constant Cp that is appropriate for each period range 
P (i.e. EHJ, VHJ or HJ). The model distribution is then given by 

"^^^ CpM'^. (16) 



dlnM 

where 



max mm 



and fp is the fraction of stars with planets in period range P with masses between Mmin = 0.1 
and Mmax = 10. 

We assume a simple planetary mass-radius relation of the form 

with 7 = 0.5, Rq = 1.5, or 7 = 0.38, Rq = 1 (see figure [I9]). We consider two different 
mass-radius relations to determine the sensitivity of our results to the assumed relation. 
The expected number of planet detections in the EHJ, VHJ and HJ period ranges is then 
given by: 

(19) 

where the sum is over all stars in the sample, Pdet,p,i{R) is the planet detection probability 
for star i for planets with radii R restricted to the period range P, and Rmin is a minimum 
planetary radius below which P^et = for all stars (we take Rmin = 0.3). Note that we 
assume that the planetary mass and period distribution is independent of stellar mass. 
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Fig. 19. — The masses and radii of solar system objects (filled circles) and transiting planets 
(open circles) are compared with two simple mass-radius relations that we adopt for estimat- 
ing the expected planet yield. The solid line is a power-law of the form R/Rj = 1.5{M/Mj)^'^ 
for M < Mj and R = 1.5Rj for M > Mj, the dotted line is a power-law of the form 
R/Rj = l.OfM/ Mj)"-^^ for M < Mj and R = R j for M > Mj. The transiting planet data 



was taken from [Torres. Winn and HolmanI (120081 ). For comparison we also show theoreti- 
cal mass-radius relations for a pure ice planet (dashed line), a gas giant with a 25Mffl core 
(long dashed line), and a pure gas giant (dot-dashed line) from iFortney. Marley. and Barnes 
( 120071 ). The gas giant planet relations are for a planetary semi-major axis of 0.045 AU. 
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We expect that the sample of stars surveyed by the RV planet searches, which is bi- 
ased towards high metallicity, provides a better match to the metallicity of the cluster than 
the stars surveyed by the OGLE transit search, whereas the OGLE sample provides a bet- 
ter match to the metallicity of the field stars than the RV sam ple. We therefore use the 
frequencies of HJ and VHJ plane ts from ICumming et al.l (120081 ) for the cluster and the fre- 
quencies from I Gould et al.l (120061 ) for the field. We take the frequency of EH J inferred from 
the SWEEPS survey for both the field and the cluster. Fixing Chj, Cvhj and Cehj for 
cluster members such that the fraction of stars with HJ, VHJ and EH J larger than O.IMq 
is 0.012, 0.0017, and 0.0005 respectively, we find Chj = 2.4 x IQ-'^, Cvhj = 3.4 x lO"^ and 
Cehj = 1-0 x 10~^. For the field stars we use the frequencies from iGould et al.l (120061 ) to 
set Chj = 6.2 x 10"^ and Cvhj = 2.8 x 10"^ 

Using the above model and Pdet{R) values from selection criteria set 2, we find that 
for cluster members we would expect to detect ~ 0.10 — 0.12 EHJ planets, 0.11 — 0.16 
VHJ planets, and 0.23 — 0.35 HJ planets where the range depends on the assumed mass- 
radius relation. For field stars, on the other hand, the expected number of detections are 
~ 0.34 — 0.51, 0.30 — 0.54 and 0.18 — 0.39. We conclude that for the above model we would 
have expected to detect ~ 1.3 — 2 stars in our entire s urvey, and theref o re our observations 
are consistent with the model. Finally, we note that iBeatty fc Gaudil (120081 ) predict that 
there are ~ 6 transiting HJ and VHJ planets per square degree orbiting Sun-like stars with 
< 20 at Galactic latitude b = 3.1° (see figure 8 of that paper). For our 0.16 square- 
degree survey they would predict ~ 1 transiting planet detection, which is comparable to 
the expected planet yield around field stars from our simulations. 



7. Discussion 



We have presented the results of a deep ~ 20 night survey for transiting hot planets in the 
open cluster M37. This survey stands out from previous ground-based transit surveys both in 
terms of the size of the telescope used, and the photometric precision attained. We observed 
~ 1450 cluster members with masses between 0.3 Mq ^ M < 1.3 Mq as well as 7814 Galactic 
field stars with masses between 0.1 Mq < M < 2.1 Mq. While no candidate planets were 
found among the cluster members, we did identify one candidate extremely hot Jupiter with 
a period of 0.77 days transiting a Galactic field star. However, the follow-up spectroscopic 
observations needed to confirm the planetary nature of this candidate would be difficult 
(although perhaps not impossible) to obtain with current technology, given the faintness 
of the source. We note that if this candidate is real, then we conclude that 0.09]'^o;q77% of 
FGKM stars have Jupiter-sized planets with periods between 0.5 and 1.0 days. This result 
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would be consistent with the results from the SWEEPS survey, and would confirm this new 
class of ultra short period planets. Wc also note that this planet frequency is small enough 
that these planets could have escaped detection in most other radial velocity and transit 
surveys. 

The primary result of this survey is an upper limit on the frequency of planets smaller 
than l.ORj. For cluster members, wc find that at 95% confidence < 25% of stars have planets 
with radii as small as 0.35-Rj and periods shorter than 1.0 day, < 44% of stars have planets 
with radii as small as 0A5Rj and periods between 1.0 and 3.0 days, and < 49% of stars have 
planets with radii as small as 0.6Rj and periods between 3.0 and 5.0 days. The upper limits 
on the smallest planets may be as much as a factor of 50% higher if all the variable stars 
near the cluster main sequence are cluster members. For the field stars we are able to place 
95% confidence upper hmits of 16% on the fraction of stars with planets as small as O.SRj 
with periods less than 1.0 day, 8.8% on the fraction of stars with planets as small as 0.5Rj 
with periods between 1.0 and 3.0 days and 47% on the fraction of stars with planets as small 
as O.bRj with periods between 3.0 and 5.0 days. Wc estimate that these upper limits may 
be higher by at most a factor of ~ 11% due to binarity. While these limits do not approach 
the observed frequency of Jupiter-sized planets with similar periods, they do represent the 
first limits on the frequency of planets as small as Neptune. We can now state empirically 
that extremely hot Neptunes (periods shorter than 1 day) are not ubiquitous, nor are very 
hot planets with radii intermediate between Neptune and Saturn. 

The limits that wc place on Jupitcr-sizcd planets arc more stringent than previous open 
cluster transit surveys, but are still above the frequencies measured by RV and Galactic field 
transit surveys. The primary limitation on open cluster transit surveys appears to be the 
paucity of stars in these systems. To place a hmit on the frequency of HJ that is less than 
2% with the same set of observations, M37 would have to have been ~ 4 times richer than 
it is. We also note that for a relatively young cluster like M37, variability may reduce the 
detectability of Neptune-sized planets by as much as ~ 50%. 

We are grateful to C. Alcock for providing partial support for this project through 
his NSF grant (AST-0501681). Funding for M. Holman came from NASA Origins grant 
NNG06GH69G. We would like to thank G. Fiiresz and A. Szentgyorgyi for help in preparing 
the Hcctochcllc observations, C. Burke for a helpful discussion, and the staff of the MMT, 
without whom this work would not have been possible. We would also like to thank the 
anonymous referee for several suggestions which improved the quality of this paper, and the 
MMT TAG for awarding us a significant amount of telescope time for this project. This 
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Table 1. Parameters from fitting tlie noise model in equation [2] to the median RMS-r 
relations after binning the light curves on several time-scales. 



Binning Timo-Scalo (minutes) 
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''The r magnitude of a source which has equal red noise and white noise 
when the light curve is binned at the specified time-scale. 



Table 2. Candidate Transiting Planets 
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20, 


.64 


23, 


.78 


22.26 


20. 


17 


111000111 


2 


2.027186 


53726.18244 


0, 


.1789 


0, 


.025 


26.33 


170049 


05: 


:52: 


: 44.65 


+32:25:55, 


.3 


18, 


.36 


17, 


.65 


17, 


.30 


18, 


.74 


17.90 


16. 


,87 


000000001 


2 


2.261081 


53725.22479 


0, 


.0181 


0, 


.055 


10.29 


170100 


05: 


:52: 


:38.44 


+32:23:29, 


.5 


20, 


.34 


18, 


.73 


17, 


.80 


21, 


.02 


19.43 


17. 


,34 


000001000 


3 


2.746378 


53725.19817 


0, 


.2963 


0, 


.025 


20.78 


230082 


05: 


:52: 


:01.32 


+32:32:51, 


,0 


17, 


.71 


17, 


.08 


16, 


.71 


18, 


.15 


17.36 


16. 


,27 


011000011 


2 


3.335057 


53726.43670 


0, 


.0280 


0, 


.025 


15.05 


270149 


05: 


:51: 


:59.04 


+32:20:38, 


.0 


20, 


.59 


19, 


.47 


18, 


.87 


21, 


.68 


19.98 


18. 


,42 


000000111 


3 


1.097491 


53724.89160 


0, 


.2132 


0, 


.065 


11.11 


280287 


05: 


:51: 


:37.41 


+32:45:03, 


.0 


23, 


.80 


22, 


.12 


20, 


.78 


24, 


.68 


22.82 


20. 


,29 


111000011 


2 


0.456034 


53725.41220 


0, 


.1331 


0, 


.055 


17.35 


330224 


05: 


:51: 


;30.37 


+32:30:29, 


,3 


21, 


.64 


20, 


.46 


19, 


.91 


22, 


.11 


20.94 


19. 


.46 


010000110 


1 


0.648749 


53725.20636 


0, 


.0502 


0, 


.060 


13.79 



CO 



''The magnitude of the source in the photometric catalog (Paper I) 
''Value from lKalirai et al.l bOOll l 

'^Value from transforming ri photometry to the Iq photometry of iNilakshi fc Sagad l l2002h 

'*9-bit flag to indicate which selection criteria sets and resolutions selected the candidate. From the left, the first three bits indicate if the candidate was selected 
by criteria sets 1, 2 and 3 respectively, using a high resolution BLS search. The middle three bits are for the low resolution BLS search applied to candidate cluster 
members, while the last three bits are for the low resolution BLS search applied to field stars. 

°This integer indicates if the BLS parameters in this table were determined for: 1. the fully processed, trend-filtered light curve, 2. the light curve processed 
through steps 1 and 2 of ^ 13. II only, or 3. the un-processed light curve. 

'Signal to Pink noise ratio (eq. |4]l. 



Table 3. 95% Upper limits on the Planet Occurrence Frequency 



Period 


Selection 


0.3Rj 


0.35Rj 


OARj 


0A5Rj 


0.5ilj 


O.GRj 


0.7ilj 


O.SRj 


0.9iij 


l.ORj 


1.5Rj 


2.0iij 


Oluster JiHJ 


1 


1.000 


0.421 


0.139 


0.067 


0.042 


0.023 


0.017 


0.015 


0.014 


0.014 


0.026 


0.048 




2 


1.000 


0.247 


0.091 


0.049 


0.032 


0.019 


0.015 


0.012 


0.011 


0.011 


0.012 


0.025 




3 




















0.014 


0.010 


0.014 


■\7XJ T 
V-tlJ 


1 


1.000 


1.000 


1.000 


0.754 


0.328 


0.117 


0.067 


0.048 


0.040 


0.036 


0.050 


0.100 




2 


1.000 


1.000 


1.000 


0.444 


0.213 


0.086 


0.051 


0.037 


0.031 


0.027 


0.023 


0.056 




3 




















0.044 


0.026 


0.031 


HJ 


1 


1.000 


1.000 


1.000 


1.000 


1.000 


0.967 


0.416 


0.248 


0.186 


0.153 


0.165 


0.300 




2 


1.000 


1.000 


1.000 


1.000 


1.000 


0.485 


0.226 


0.139 


0.104 


0.083 


0.069 


0.149 




3 




















0.202 


0.104 


0.124 


Field EHJ 


1 


0.230 








0.023 




0.007 






0.004 


0.003 


0.004 




2 


0.158 








0.014 




0.005 






0.003 


0.002 


0.002 




3 




















0.003 


0.002 


0.002 


VHJ 


1 


1.000 








0.148 




0.035 






0.012 


0.007 


0.007 




2 


1.000 








0.088 




0.020 






0.008 


0.005 


0.004 




3 




















0.009 


0.005 


0.005 


HJ 


1 


1.000 








1.000 




0.240 






0.062 


0.028 


0.026 




2 


1.000 








0.469 




0.095 






0.027 


0.013 


0.012 




3 




















0.044 


0.020 


0.019 
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Table 4. 95% Upper limits on the Planet Occurrence Frequency for Field Stars Using 

Different Galactic Models 



Radius [Rj] EHJ'' VHJ HJ 



Trilegal Model 





.3 


0. 


,143 












0, 


.5 


0. 


.014 


0, 


,082 


0, 


,401 







.7 


0. 


,005 


0, 


,022 


0, 


,104 




1, 


.0 


0, 


,003 


0, 


,009 


0, 


,033 




1, 


.5 


0, 


,002 


0, 


,005 


0, 


,015 




2, 


.0 


0, 


,002 


0, 


,005 


0, 


,013 


Besangon Model 


0, 


.3 


0. 


,151 










Ay = 0.5 mag/kpc 


0, 


.5 


0. 


,010 


0, 


,065 


0, 


,394 




0, 


.7 


0. 


,004 


0, 


,016 


0, 


,070 




1, 


.0 


0. 


,003 


0, 


,007 


0, 


,023 




1, 


.5 


0, 


,002 


0, 


,005 


0, 


,014 




2, 


.0 


0. 


,002 


0, 


,005 


0, 


,015 



^Tlic 95% upper limit on the occurence frequency of planets 
with the specified Radius and within the EHJ period range. 



Table 5. 95% Upper limits on the Planet Occurrence Frequency for Cluster Stars 
Assuming All Variables are Cluster Members 



Radius [Rj] 


EHJ^ 


VH.I 


HJ 


0.35 


0.381 


1.000 


1.000 


0.40 


0.123 


1.000 


1.000 


0.45 


0.060 


0.633 


1.000 


0.50 


0.037 


0.281 


1.000 


0.60 


0.021 


0.103 


0.636 


0.70 


0.015 


0.058 


0.277 


0.80 


0.013 


0.040 


0.161 


1.00 


0.011 


0.028 


0.090 


1.50 


0.012 


0.024 


0.071 


2.00 


0.022 


0.056 


0.154 



"The 95% upper limit on the occurence 
frequency of planets with the specified Ra- 
dius and within the EHJ period range. 



